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Errata  for  AM^  Note  No.  16 

Two  minor  computational  ©rrora  which  do  not  affect 
the  validity  of  the  report  have  been  found  in  this 
Note,  as  follows: 

Page  2,  third  line  from  bottom: 

Instead  of:   7.581x^  +  4.053x®  +  x^ 

Read:        7.481x^  +  4.0553c^  +  x* 

Page  10,  fourth  line  from  top: 

Instead  of:  x  +  ,725x  +  1.98 
Reads        x^  +  .725x  +  1.99 
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INTRODUCTION 

V/liile  the  present  note  doe:;;  not  originate 
from  a  specific  request,  its  publication  by  the 
Applied  Mathematics  Panel  seems  justified  since 
it  concerns  the  problem  of  determining  numerical- 
ly the  complex  roots  of  algebraic  equations.   With- 
in the  NDRC  and  in  rslated  British  agencies,  the 
need  for  such  numerical  iriethods  has  become  apparent^ 
hence  the  follov/ing  memorandum  by  Drt  B.  Friedman 
of  the  Nov;  York  University  Applied  Mathematics  Group 
is  being  issued  by  the  Applied  Mathematics  Panel. 
This  note  is  an  amplification  of  a  short  abstract 
presented  to  the  American  Mathematical  Society 
(Bulletin  of  American  Mathematical  Society,  v.  49, 
page  859,  1943). 

R.  Courant 

Contractor's  Technical 
Representative 
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COL'iPLEX  ROOTS  OF  POLYNOMIAL  EQUATIOIJS 

The  problem  of  finding  the  complex  zeros  of  a 
polynomial  equation  is  important  in  many  branches  of 
science  and  engineerinn;.   Graef  f e '  s  method,  which  is 
very  powerful  when  the  roots  are  I'ual,  is  not  practical 
when  most  of  the  roots  are  complex. 

Recently  a  British  report"  described  a  method 
for  calculating  the  roots  of  equations  of  the  fifth  and 
sixth  degree.   One  disadvantage  of  the  proposed  method 
is  that,  at  each  stop,  it  requires  the  solution  of  a 
quadratic  equation. 

This  report  presents  a  simple,  general  method  for 
finding  the  complex  roots  of  any  polynomial  equation. 
All  the  operations  can  be  easily  performed  on  a  calculating 
machine. 

Since  complex  numbers  are  awkward  to  handle,  it  is 
desirable  to  operate  vi/ith  the  quadratic  factors  of  the 
polynomial  instead  of  ;.vith  its  complex  roots.   The  question 
of  finding  quadratic  factors  can  be  divided  into  two  pror 
blems:   First,  to  obtain  quickly  an  approximate  value  for  the 
quadratic  factor  and  second,  to  inprovo  the  approximation. 
The  first  problem  will  be  answered  by  a  method  of  successive 


■:;-  V/A-1907-1  A.D.R.D.E,  Res.  Rep.  Wo.  238 
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forward  and  back  divisionsj  the  second,  by  an  extension  of 
Nev/ton's  method  to  complex  roots. 

Method  I«  Succeasive  Division 

The  method  can  be  beat  described  by  considering  a 
par ti cular  example : 

Find  the  roots  of 

f(x)    =  x^   +  5x^  +  13x^   +  15x  +   10  =    (x^   +  2x   ^  2){x^'   +  3x   +  5) 

2 

Assmne   as   a   trial  factor     x      +  Ox   +   .25    .    (A  method 

for   determining  the   first   trial   factor  will  bo   given  later. ) 
Then 

/(^^        =  x2   ^  5x   +   12,75  ^   H«'^5x   ^  6.8125   ^ 
x^   +   .25  x^   +   .25 

2 
Divide      f(x)   by   12.75   +  5::   +  x      in  ascending  powers    of     x. 


as  follows: 


2 


.784   +   .947x   -f   .595x^ 


12.75   +   5x   +  X  r~~~7  ,    T  „   2  .    r-    3  4 

'     10  +   16x  +  ISx'^  +  5x        +  X 

10   +      5.92X   +      .784x^ 

12.08X   +12.216x2  +   5x^ 

12.08X   +  4.735x^  +      .947x^ 

7.581x2  -f  4.053x'^   +  x^ 

7.581x2  +  2.975x^   4-   .595x^ 

l.lSOx^   +   .405x^ 


-"-  This    is   an  extension  of   a  method  proposed  by  Lid,    Journel 
of  Mathematics   and     Physics,    June    1943. 
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2 
The   quotient   is    .SQSx     +   .947x  +   .784   or 

.595    (x^   +   1. 59x   +   1.32),    so   that     x^   +   1.59x   +   1.32      is 

2 
a   closor  approximation  to   the   factor     x     +  2x  +  2. 

The   motivation  for    this   method  can  bo   scon  from 

the   following  facts: 

p 

Suppose  the  correct  factor  x  +  2x  +  2  v.'cre  knov;n. 

Then 

~^^^ =   x2  ^  3x  ^  5 

x*^  +  2x  +  2 

and  if  f(x)  wore  divided  by  5  +  3x  +  x    (in  ascending 
powers  of  x),  the  quotient  would  be  x'^  ■♦■  2x  +  2,    the  correct 
factor.   This  suggests  that  an  iteration  method,  which  con- 
sists in  dividing  f(x)   by  trial  factors  in  descending 
and  ascending  powers  of  x,  successively,  may  converge  to 
the  correct  factor. 

The  method  does  converge  for  the  example  given  pre- 
viously.  After  five  repetitions  of  the  process,  the  result 

2 

is   X  +  1.98x  +  1.99.   By  further  repetitions,  any  desired 

accuracy  can  be  obtained. 

The  work  necessary  for  this  method  can  be  simplified 
if  synthetic  division  is  used.   Suppose 

f(x)  =  x'^  +  a^x^^*"^  +  agX^"^  +  , . .  +  a^ 

and  the  trial  factor  is   x""  +  px  +  q.   Use  only  the 
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coefficients  and  perform  the   ivision  as  follows: 


1 
-P 


+  a,   +  ap   +  a^  +  a^   +  •..  +  a 

-   q    -  b^q  -  bgq  -  . . ,  -,  b^^^q 


1  +  b^  +  bg  +  b^  +  b^   +  '••  +  ('in'"^n-2^^ 


whcr( 


b,  =  a^  -  p 


b^  =  a^  - 


b/,  =  a„  - 


'4 


b;^P 
bgp 
b3P 


q 

b^q 
bgq  etc* 


Notice  that  the  signs  of  the  coefficients  of  the 
quotient  have  been  changed*   This  permits  us  to  add  instead 
of  subtracting.   The  numbers  in  the  successive  diagonals 
are  obtained  by  multiplying  the  coefficients  of  the  quotient 
by  the  coefficients  of  the  divisor* 

As  further  illustration,  suppose  we  try  the  ne:;:t  step 
in  this  method  for  the  example  mentioned  previously,  that  is, 
divide  f  (;■:)  by  x^   +   1.59::  +  1«32, 

1 

-  1.59 

-  1.32 

1  +  3.41  +  5,258  +  1.55  +  1.14 


1 

+  5       +13        + 

16       + 

10 

-1.59      -   5.422 

^9.95 

-   1.32 

-4.50 

-8.26 
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The   quotient   Is     x"^   +  3.4lx    +  6.258     while    the 
remainder   is      l,55x    +  1.14-.      The   next   step  would   be    to 
divide      10   +   16x    +  13x^   +   5x^   +  x'*      by      5.253    +  5.4Lx    +  x^. 
The   work   is   simplified    if    the    leading    coefficient   of   the 
divisor   is    one.      V/rite      S,253    +  3.41j:    +  x        as      6.258  times 
(1   +   .54x    +  ,16x    ).      The   division   proceeds    as    follows ; 


.16 


10    +    16 

+   13 

+   5    +    1 

-      5.4 

-  5.7 

-  l.G 

10    +   10,6      +      5.7 


Since    the   remainder   is    not   needed,    the   work  can   stop  here* 

How   the   quotient   is      10   +  10-,6x    +   5.7x"^   =   5.7(x'^   +   l,86x   +  1,75) 

2 

and  X   +  1.86x  +  1.75   is  the  approximation  to  the  correct 

factor. 

The   whole    proce.cs    can   be    combined   as   follows;''*' 


1 
1.59 

1,32 


1 

+ 

5 

+ 

13 

- 

1. 

59 

- 

5, 
1. 

422 

32 

1   +  3.41   +      6.258 


13         +15        +   10 

-  5,7    -      5»4 

-  1.6 

+   5.7    +    10.6    +    10 


-  i54 

-  .16 


Kote    that    6,258(1   +  ^''^L-"    +  ^■■o^a'  ^^)    =  6.258(1 +,54x  +  ilSx^ ) 


5.^53- 


The   next   approxim.ation.  is  x      +      h'w  x   +  ~^,  „   . 

Two  more    repetitions    of   this   forward    and   back  di- 

2 
vision   give    as    the   factor     x      +  l,98x    +  1.99,      The    process 

can  be   continued   to  give   any  needed   accuracy.. 


The   preceding   synthetic   division  is   written  in  the 
reverse   direction  merely  for   convenience   of   computation, 
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This  method  can  be  applied  to  polynomials  of  any 
degree.   In  general  it  v/ill  converge  to  the  quadratic 
factor  containing  the  two  roots  of  smallest  absolute 
value.   The  initial  trial  factor  should  be  chosen  so  that 
its  roots  arc  smaller  in  absolute  value  than  those  of  the 
given  equation. 

.  Simple  reasoning  gives  a  lower  bound  for  the  ab- 
solute value  of  the  roots.   Consider  the  equation  gi.ven 
in  the  British  report  referred  to  on  page  1: 
g(x)  =  X®  +  14.71x^  +  92.76x^  +  273. 3x"^  +  5l5x^  +  532z 

+  435.5  =  0  . 
The  following  inequality  is  obtained: 

|x^  +  I4.71x^  +  92.7Gx'^  +  273. 3x"^  +  515x^  +  532x  +  435.51 
>  435.5  -  532|xl  -  515|x^I  -  273.3|x^|  -  92.76[x'^I 
-  14.7l|x^|  -  IxT'^ 

The  right  hand  side  is  greater  than  zero  if  |x|  >  .77. 

2        2 

A  trial  factor  can  be  then  taken  as  x  +  (.77)   or 

2  .   . 


The  work  proceeds   as   follows; 


1 

+ 

14, 
0 

71 

+ 

92. 
0 

76 

+ 

273. 

0 

3 

+ 

515 
0 

- 

. 

G 

- 

8. 

8 

- 

55. 

3 

515 

+ 

532 

+ 

435, 

5 

-162 

- 

248 

-   87 

1 

-.57 
-.2 


1   +   14.71   +   92.16    +   264.5   +  459.7 


Note    that   the    "back"   divisor   is    1  +  -1^4^  :■'   '^A^^'i^  x^ 


=   1   +   •57x  +   ,2: 


266    +   234   +   435.5 

64.5        

459.7  '■■    '459.7 
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The    second  approximation  is     x^   +  |§§x+  ^g^'P 


X     +   l.Lx   +   1.6 


1 

1   +    14.71   +   92.76    +   273.3   +   515 

-1.1 

-      1.1      -    14.97    -      83.8  -    134.4 

-1.6 

-      1.6      -      21.3   -    121.9 

515        +   532        +  4  35.5 
-146.9   -    348.4 
-161.1 


1 

-  .80 

-  .37 


1   +   13.51   +  76.2      +    167.7   +   208.7  207        +   183.6    +  455.5 

The    third  approxination  is     x^   +  l^.S'^"^''   +  ^^^'^ 


207 


207 


=  X      +   .89x   +   2.1    . 


1   +    14.71   +   92.76   +  273.3   +   515  515        +   532        +  435.5 

.89    -    12.30   -      69.7    -   155.4         -126         -   392 
-      2.1      -      29.0   -   164.6         -174.2 


1 

-  .90 

-  .40 


1   +    13.82   +  78.36    +   174.6   +   195.0  214.8   +   140        +  435.5 

2 

The   fourth  approximation  is      x      +   •G5x   +   2.03   . 

This  method   can  be   continued  as    long  as   necessary, 
but,    since   the  approximation   is    fairl;;-   close    to   the    correct 
factor,    another  method  which  converges   more   quickly  can  bo 
used. 

Method   II.    Generalized   rJcwton's    IvlGthod" 

Suppose     f(x)      is   divided  by  some    trial   factor 
X     +  px   +   q.      A  quotient      Q(x)      and  a  remainder     rx   +   s 
v;ill  be   obtained: 


A   somewhat   similr.r  method   has   been  proposed   independently 
by  Hitchcock,    Journal  of   Mathematics   and  Physics, 
May   1944. 
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2 
(1)      f(;;)    =    (x      +  px   +  q)      Q(x)    +  r;-:   +   s,    where      r     and 

s   are  functions  of  p   and  q.   The  problem  is  to  change 
p   and  q   so  that  r   and  s   will  be  zero.   Using  Tay- 
lor's theorerii  in  two  variables,  this  can  be  done  if  the 
partial  derivatives  of   r  and  s  with  respect  to  p   and 
q  are  Imown. 

Differentiate  equation  (l)  with  respect  to   p.   The 
result  is 

(3)   0  =  .xq(x)  .  {:/  ^  px  +  q)  |S  ^  ||x  +  3| 

2 

if  the  zeros  of  x   +  px  +  q   are  substituted  for  x, 

it  is  seen  that 

^x  +  ^  =  remainder  of  [-xQ(:c)]  when  divided  by 

2 

X   +  px  +  q. 

Similarly 

^x  +  ^  =  remainder  of  [-Q(x)]   when  divided  by 

2 

X   +  px  +  q. 


Using  Taylor's  theorem,  the  solution  of  the  equations' 
(3) 


||Ap.||Aq.-r 


op      aq   ^ 
will  give  the  changes  in  p   and  q  which  will  cause  r 

and  s   to  become  zero. 

As  an  illustration  consider  the  previous  example. 

f(-:)  =  x*^  +  I4.71x^  +  92.76x'^'  +  273. 3x^  +  515x^ 

+  532x  +  4  35.5 

2  P 

where  x   +  px  +  q  =  x  +  .65x  +  2.03 


for  A  p   and   Zi^q. 
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1 

.   .65 
•2.03 


1 

-   .65 
■2.03 


1 

■   .65 
-2.03 


(4) 
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1   +    14.71   +   92.76 

+   273.3   +   515        +   532        +  435.5 

.65    -      9.14 

-      53.0   -   124.7    -   146. 1 

-      2.03 

-      20. 5    -   165.6    -   389.4    -   456. 1 

1   +   14.06    -f   81,59   +   191.8   +   224.7    -        3.5    -      20.6 
rv   +   s   =   -3.5;j   -  20.6 


■^"   +    14.06::^   +   81.59x^   +    191. 8x   +   224.7 


1 

+ 

14 

.06 

+ 

81. 

59 

+ 

191. 

.8 

+ 

224. 

7 

-- 

.65 

- 

8. 

72 

- 

46. 

.05 

- 

2. 

03 

- 

27. 

.22 

- 

143. 

8 

1   +   13.41   +  70.84   +   118.53   +      80.9 
Therefore  ^  =   -118.53 


as 


-   80,9 


To   find  the  partial  derivatives  with  respect   to  p, 

2 

divide  ycQ(x)      by     x     +   .55;-:  +  2.03.      Note    that  most    of   the 

work   is    identical  v/ith   that  in   the  preceding  division. 


1 

+ 

14, 

,06 

+ 

81. 

59 

+ 

191. 

8 

+ 

224, 

7 

+ 

0 

- 

. 

65 

- 

8. 

.72 

- 

46. 

05 

- 

77. 

0 

- 

2. 

,03 

- 

27. 

22 

- 

143, 

_8_ 

_- 

240. 

_5 

1   +   13.41   +  70,84   +   118.53   +        3. 9    -   240.5 

Then||  =  -3.9    ,      ||.     =  +  240.6 
Equations    (3)   become 
-   3.9Ap    -   118.53Aq  =   3.5 
240. 6  Ap    -      80.9  Aq   =   +20.5 
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Their    solution   is 

^p   =   +    .075    ,     Aq   =    -    .032 
so   that   the   next  approximation  to   the   quadratic  factor    is 
x^   +   .725x   +   1.98. 

One  more   stop: 


■   .725 
•1.98 


1 

a. 

14 

.71  + 

92. 

76 

+ 

273. 

3 

+ 

515 

+ 

532 

+ 

435. 

5 

- 

.725- 

10. 

14 

- 

58. 

46 

- 

135. 

68 

'•* 

159. 

25 

- 

1. 

_98 

- 

27. 

_69_ 

- 

J5_9. 

67 

- 

370. 

56 

- 

434. 

91 

1  +    13.985-   80.64    +    187.15   +   219.65   +        2. 19   +        0.59 


r   =   2.19 


s    =   0.59 


1 
•    .725 

-1.98 


1   +    13.985   +   80.64   +    187.15   +   219.65   +   [O] 

.725    -      9.61   -      50.06    -   [80.36] 
-      1.98   -      26.25    -    156.72    -    [219.46] 


1   +   13.26      +   69.05   +   110.84    +      82.93    -    [219.46]    "' 

+   [    2.57] 

Equations    (3)   become: 
-   2.57Ap    -   110.84Aq   =   -  2.19 
219.46Ap    -      82.93Aq   =   -   0.59 
Their    solution   is 
Ap    =    ,0047    ,         Aq   =    .0195 


and  therefore   the   factor    is      x' 


.7297X   +   1.996.      The   othc] 


factor    is      x'^   +   13.98x^  +  80.56x^  +   186. 6x   +  218.      This    can 
be   factored   in  the    same   way  as    the   original  pol7/nomial. 

Add!  t  i  o  na  1  T'lo  t  e  s 

1.   If  all  the  roots  of  the  polynomial  equation  have 
their  real  parts  negative  (this  is  the  most  coim-.ion  case  in 
practice),  then  a  better  first  approximation  can  be  found 


In  order  to  combine  the  division  of  Q(x)  and  xQ(x)  into 
one  operation,  the  additional  nuinbers  needed  are  put  into 
brackets.  RESTRICTED 
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by  using  an  idea   of   Lin  in   the  paper   previously  mentioned. 
Let 
f(x)    =  x""  +  a^x""--   4-  ^2^^"^   +    ...    +   a^^_^x   +  a^ 

Use   as   a   first  approximation 

^  I    ^   ^  ^n-1  ..   ^  ^n      , 

a      rjX      +  a.^     X   +  a      =  a      pix      + ^   + ; 

n-2  n-1  n  n-<d  n-2  n-2 

2.   Method  I  will  converge'  to  the  quadratic  factor 
whose  roots  are  sniallest  in  absolute  value.   For  an  equa- 
tion of  odd  degree  vvith  one  real  root  smaller  in  absolute 
value  than  all  the  other  complex  roots,  it  is  necessary  to 
remove  the  real  root  before  using  Method  I. 

Another  possibility  is  to  find  the  quadratic  factor 
whose  roots  are  largest  in  absolute  value.   This  can  be  done 
by  assuming  a  trial  factor  and  then  dividing  in  back  first, 
that  is,  divide  in  ascending  pov./ers  and  then  in  descending 
powers,  instead  of  dividing  in  descending  powers  and  then 

in  ascending. 

3.   At  each  step  in  Method  I  the  difference  between 
the  coefficients  of  the  exact  factor  and  the  coefficients 
of  the  trial  factor  is  approximately  k  times  the  corres- 
ponding difference  at  the  preceding  step.   k  depends  . 
directly  on  the  ratio  of  the  absolute  value  of  the  smallest 
pair  of  roots  to  the  absolute  value  of  the  other  roots. 

At  each  step  in  liethod  II  the  coefficient  difference 
is  approxime.tely  the  square  of  the  corresponding  difference 


A  proof  of  convergents  will  be  given  in  a  mathematical 
journal  at  some  future  date. 
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1 

-2.9 
-4.2 
■2.4 


at  the  preceding  step.   Therefore,  if  the  first  approxina- 
tion  is  close  enough.  Method  II  will  converge  :nuch  more 
quickly  than  I. 

In  practice,  if  more  than  5  per  cent  accuracy  Is 
needed  in  the  roots.  Method  II  should  bo  used  after  four 
or  five  trials  of  Method  I. 

4.   Both  methods  can  be  applied  to  find  real  fac- 
tors of  any  degree.   As  an  example  we  shall  find  a  cubic 
factor  of  x^   +  IVx'^'  +  lOSx"^  +   297x^  +  422x  +  240. 

Per  a  first  approximation  take 

103x^  +  297x   +  422x  +  240  or  x   +  -^  x   +  ^  x  +  -^  . 

lOo  IQo  lOo 

The   first   step    in  Method   I   can  be   sot  up  as   follows t 


1+17         +    103  103        +   297        +  422        +  240 

-      2.9    -      40.9         -  49.7    -      88.4    -      53.3 
4.2        -      6.2    -        4.1 


1 

-  .243 

-  .017 


1   +    14.1   +      57.9 


47.1   +   204.5   +   363.7   +   240 


The   second  approximation   is 
.„    T    /    3    ,    204.5      2    ^    363.7 

47.1     (X        +    -  yi,-,    ■-,     X        + 


T7. 1  '^ 


'vm 


X    + 


MO 
"5771 


)    . 


After   five    steps    the  result   is 


(x^   +   5.5x'^   +    11. 3x   +   8.3)    (x      +    11. 5x   +   23.5)    . 
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